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I. INTRODUCTION 


The Goddard Laboratory for Atmospheres (GLA) has over the past few years 
acquired a large data base of 00Z and 12Z initialized analyses (1980-87) from the European 
Center for Medium Range Weather Forecasts (ECMWF). In support of our on-going 
research activities we have found it convenient to compile these analyses in the form of an 
atlas. 

The original ECMWF/WMO data was obtained on a 2.5 ° latitude by 2.5 0 longitude 
grid in a packed format. This data was unpacked and interpolated to a 4 0 latitude by 5 ° 
longitude grid for the purpose of computing various statistics and plotting. Appendix A 
describes further details of the data and data processing. In part I of this atlas we present 
various first moment quantities for monthly, seasonal and annual averages. The present 
volume (Part II) presents various second moment quantities averaged over a season 
including both transient and standing eddy components. Since we accumulate the statistics 
on a monthly basis, the majority of the transient statistics presented in this volume are 
computed as deviations from monthly means. However, selected fields have been 
processed differently in order to retain the longer time scales and remove the influence of 
the seasonal cycle. These statistics are based on deviations from the average seasonal 
cycle. The deviations are filtered to separate the contributions of the band pass (2.5 - 6 
days) and low pass (> 10 days) regions of the spectrum. 

The ECMWF analysis system has undergone a number of changes during the 
period 1980-87. Trenberth and Olson (1988) have done an extensive study of the ECMWF 
analyses for 1979-86 and examined the impact of analysis system changes on data quality 
(see their Table I and discussion). The change to diabatic nonlinear normal mode 
initialization in September of 1982 dramatically increased the strength of the tropical 
divergence. For quantities sensitive to the divergence we have limited the time averages to 
1984-87. For these fields the shorter time-averaging period is clearly indicated. Otherwise 
the averages include all 8 years. It must be emphasized that many other changes occurred 
to the analysis system during this time period leading to inhomogeneities in the data. This 
is particularly true for the moisture fields (see eg. Trenberth and Olson, 1988). No attempt 
has been made to take this into account when performing the time averages. 

We note that other recent compilations of assimilated data that may be useful for 
comparisons are those of Trenberth and Olson (1988) based on essentially the same data 
(1979-86), Hoskins et. al. (1989) which emphasizes special diagnostics and year to year 
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variations in ECMWF analyses for the period 1979-1989 and Lau (1984) which examines 
the FGGE year based on analyses produced at GFDL. 

In accordance with ECMWF conditions any publications making use of the figures 
and/or data in this adas must acknowledge ECMWF as the data source. 


II. LIST OF QUANTITIES 


The bulk of the statistics presented in this volume are accumulated on a monthly 
basis where daily fluctuations are computed with respect to a monthly mean . This allows 
the statistics to be easily updated as new data becomes available. For the transients the 
covariances are defined as 


xy = x y + x'y' (1) 

where the bar denotes a monthly average and the prime a deviation from the monthly 
average. For zonal averages 

py] = [x] [y] + [x * y *] + py] (2) 

where the brackets denote a zonal average and the asterisk (*) denotes a deviation from the 
zonal average. The monthly mean statisti cs are then averaged to produce the seasonal and 
annual means . The global maps consist of various transient eddy covariances (last term in 
(1)), while the zonally averaged covariances consist of both the standing and transient eddy 
terms (second to last and last terms of (2)). 

For selected statistics of the form x' y' the data was processed differently in order 
to isolate particular regions of the spectrum. This was accomplished by first removing the 
time mean and annual and semi-annual harmonics computed from the average of all 8 (or in 
some cases a subset of the 8) years of data. The deviations were then further processed to 
isolate the synoptic and low frequency fluctuations. This was done by applying a band 
pass filter which retains time scales between 2.5 and 6 days and a low pass filter which 
retains time scales of 10 days and greater. The seasonal statistics were then computed from 
deviations with respect to the app ropriate seasonal mean for each year and then averaged 
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over all years . Note that only the intraseasonal time scales are retained in the deviations. 

The filters are essentially identical to those described by Blackmon and Lau ( 1980). The 
only difference is that the filters used here are of higher order (35 point instead of 10 
point). The advantage of the higher order is that it provides a better approximation to the 
ideal amplitude response function. The disadvantage is that more points are needed so that 
the seasons are to some extent more "contaminated" by the adjacent seasons. In practice, 
litde difference was found between results from the two filters. See Schubert (1986) for a 
comparison of the amplitude response functions of the present and Blackmon and Lau 
(1980) band pass filters. 

While the low and band pass filtered fields are the primary focus of the above 
analysis, for completeness, we also show the unfiltered fields. In general, the statistics 
based on the deviations from the monthly means and those from the unfiltered deviations 
from the seasonal cycle were found to be very similar. The latter are composed of the 
intraseasonal time scales (approximately 1-90 days) while the former include, in addition to 
time scales less than one month, part of the seasonal signal inherent in the deviations from 
monthly means. 

We also include the extended E-P vector 

Eu =(IM, -3v) 0) 

introduced by Hoskins et. al. (1983) for both the unfiltered and filtered deviations from the 
seasonal cycle . The present analysis uses the slightly modified form (the factor of 1/2) of 
this quantity which is parallel to the relative group velocity for the barotropic case 
(Trenberth, 1986). 

The following lists the quantities contained in Part n. The units are given as part of the 
figure captions and contour intervals are easily seen from the plots themselves. Plots are shown 
for seasonal and annual averages. The global maps are presented at 850 and 200 mb. The cross 
sections are based on all the levels obtained from ECMWF which are 1000, 850, 700, 500, 300, 
200, and 100 mb. The fields are grouped according to season. Within each group the first set of 
statistics involve deviations from monthly means (averaged over a season) which consist of both 
global fields and zonal averages. This is followed by selected maps and cross sections of 
unfiltered and filtered quantities computed from deviations from the seasonal cycle. The final maps 
are time series of zonally averaged monthly mean statistics. We indicate below the quantities and 
(for convenience) the relative positions of the plots in the atlas for each data group. 


IX 


Quantity 


VuV(200mb), VvV (200mb) 

1 

Vrr (850mb), Vt'T' (200mb) 

2 

(850mb), VZ’Z’ (500mb) 

3 

7r (850mb), uT (850mb) 

4 

Tq 1 (850mb), iTq 1 (850mb) 

5 

uV (200mb) 

6 

fflT (850mb), oTcf (850mb) 

7 

v’Z’ (200mb), u'Z’ (200mb) 

8 

^/pv], 

^f[u*u*] 

9 

JH. 

^[v*V*] 

10 

i TT i • 


11 

i z ' z i • 

^[z*z*] 

12 

V[7q’]. 

^[q*q*] 

13 


[u * v *] 

14 

[v'T], 

[v*T*] 

15 


[v*q*] 

16 

[v'Z']. 

[ v * Z *] 

17 

[0)'T], 

* 61 
* 

HI 

* r 

18 

[to'q’] , 

[w * q *] 

19 



Quantity J 

Relative Position 

Unfiltered 


v'T* (850mb), v’q' (850mb) 

20 

uV (200mb) 

21 

V Z'Z’(500mb), u’u'+v'v’ (200mb) 

22 

V'T' (850mb), V'q’ (850mb) 

23 

E u 200mb (contours of v'T' at 850mb) 

24 

[Sv], [7T]> 

25 

Low pass (periods > 10 days ) 


v'T' (850mb), v’q' (850mb) 

26 

uV (200mb) 

27 

VZ'Z' (500mb), u'u'+v’v' (200mb) 

28 

V’T’ (850mb), V'q’ (850mb) 

29 

E u 200mb (contours of v'T' at 850mb) 

30 

[uV], [v'T'] 1 ’2 

31 

Band pass (2.5 < periods < 6 days) 


v'T' (850mb), v’q’ (850mb) 

32 

uW (200mb) 

33 

VZ'Z' (500mb), u'u’+v'v' (200mb) 

34 

YT (850mb), V’q’ (850mb) 

35 

E u 200mb (contours of v' T’ at 850 mb) 

36 

[uV], [v’T]i.2 

37 


*Due to tape problems the heat fluxes (for zonal averages only) are based on the years 1980-85 
2 The zonally-averaged fluxes were processed using the Blackmon and Lau (1980) versions of the band and 
low pass filters. All other quantities are computed with the higher order versions of these filters described 
in section II 
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Quantity 


'[u’u 1 ] (200 mb) 
'[u*u*] (200mb) 
[vV] (200mb) 


v*v*] (200mb) 

[tt] (200mb) 

[T*T*1 (200mb) 


T’T (850mb) 


^[T*T*] (850mb) 

8 

^[q'q’J (850mb) 

9 

-^[q * q *] (850mb) 

10 

^[Z’Z’J (500mb) 

17 

^[Z*Z*] (500mb) 

18 

[uV] (200mb) 

19 

[u*v*] (200mb) 

20 

[v'T 1 ] (850mb) 

21 

[ v * T *] (850mb) 

22 

[Tcf] (850mb) 

23 

[v*q*] (850mb) 

24 

[v’Z] (200mb) 

25 

[v * Z *] (200mb) 

26 


Appendix B describes how to access these data sets on the IBM 3081. 
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V. APPENDIX A- THE DATA PROCESSING 


The ECMWF analyses are the initialized fields provided to GLA (code 61 1) in an 
internal bit format. The data processing involves several steps including unpacking and 
converting to a 2 X 2.5 degree grid. During this step the relative humidity is converted to 
specific humidity. Then an interpolation is made to a 4 X 5 degree grid. All eddy statistics 
are computed from the lower resolution fields. Both resolutions of the basic parameters are 
available on tape. Further details of the data are provided below. 

ORIGINAL PACKED DATA AS RECEIVED FROM ECMWF: 

2.5° latitude by 2.5° longitude in internal bit format 
1000, 850, 700, 500, 300, 200, 100 hPa 
separate 00Z and 12Z tapes 


DATA (IM,JNP) 


IM= 1, .. , 144 starting at 177.5 W and going east to 180 E 
JNP = 1, .. , 73 starting at North Pole 

(1) geopotential height (m) 

(2) temperature (°K) 

(3) u (m/s) 

(4) v (m/s) 

(5) relative humidity (%) 

(6) vertical velocity (Pa/s) 


xv 
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UNPACKED ECMWF DATA 


2° latitude by 2.5° longitude or 4° latitude by 5° longitude 
1000, 850, 700, 500, 300, 200, 100 hPa 
00Z and 12Z combined for each tape 

DATA(IM,JNP) 

IM=1, .. , 144 (72) starting at 180° W and going east to 177.5° E (175° E) 
JNP= 1, , 91 (46) starting at South Pole 

Quantity 

(1) sea level pressure (hPa) 

(2) geopotential height (m) 

(3) temperature (°K) 

(4) u (m/s) 

(5) v (m/s) 

(6) specific humidity (g/g) 

(7) vertical velocity (Pa/s) 

(8) divergence 

The specific humidity (q) is computed from the relative humidity 
q = q s * RH * 0.01 


where 


q s = 0.622 *e s /(p- .378 e s ) 

and e s is the saturation vapor pressure determined from the empirical formulae 

6.1078 * EXP { 21.8747 * (T- 273.16 )/(T- 7.66)} for T< 273.16 
6.1078 * EXP { 17.269 * ( T- 273.16 ) / (T - 35.86) } for T > 273.16 


xvi 



VI. APPENDIX B- SOURCES OF DIGITIZED DATA 


The data fields (for Part II of this Atlas) are currently stored on several MVS data 
sets and are meant for use by NASA scientists with access to the IBM 308 1 mainframe. 
We also anticipate that these fields will be made available on floppy disks for both IBM 
compatible and Macintosh personal computers in the near future. For further information 
concerning the use of this data contact 

Siegfried Schubert 
NASA/GSFC Code 61 1 
Greenbelt, MD 20771 
(301) 286-3441 

In the following we present the information necessary to access the data sets for Part II on 
the IBM 3081. 


1 . DATA SET NAMES AND READ FORMATS 
a) DEVIATIONS FROM MONTHLY MEANS 
GLOBAL FIELDS 

W3RWH.ECMWF. V ARNCE 

W3RWH.ECMWF.COVRNCE 
VT uT' v 7 ? u 7 ?’ uV 

y y > * 

W3RWH.ECMWF.COVRNCE2 
® 7 T i > ®’q , > 7Z ; ,u 7 Z 5 (1980-87) 

W3RWH.ECMWF.COVRNCE3 
®T’ > ® 7 q' ) 7z ; j u 7 Z 1 (1984-87) 

ZONALLY-AVERAGED FIELDS 

W3RWH.ECMWF.ZONMN.VARNCE 

W3RWH. ECMWF. Z ONMN . VARNCE 2 

-pz) iPH VP*’*] 

W3RWH.ECMWF.ZONMN.COVRNCE 

[uV] [u*v*] [v’T’] [v*T*] [v'q’] [v*q*] [v’Z’j [v*Z*] 

W3RWH.ECMWF.ZONMN.COVRNCE2 
[wT] [to * T *] [w 7 ^] [<o*q*] 


XV1U 



READ FORMAT 

C NSEAS IS THE NUMBER OF SEASONS = 4 
C IS=1 DJF (AVERAGED FOR 1980-87) 

C IS=2 MAM (AVERAGED FOR 1980-87) 

C IS=3 JJA (AVERAGED FOR 1980-87) 

C IS=4 SON (AVERAGED FOR 1980-87) 

C IS=5 ANNUAL (AVERAGED FOR 1980-87) 

C NLEV IS THE NUMBER OF LEVELS = 7 
C IL=1 1000MB 

C IL=2 850MB 

C IL=3 700MB 

C IL=4 500MB 

C IL=5 300MB 

C IL-6 200MB 

C IL=7 100MB 

C NVAR IS THE NUMBER OF VARIABLES IN THE PARTICULAR DATA SET 

PARAMETER(IMP1=73,JNP=46,NSEAS=5,NLEV=7,NVAR=5) 

DIMENSION VAR(IMP1,JNP,NVAR,NLEV,NSEAS) 

DIMENSION VARZ(JNP,NVAR,NLEV, NSEAS) 

C J=1 IS THE SOUTH POLE, 1=1 IS THE DATELINE 
C FIELDS ARE WRAPPED SO THAT X( 1 ,J)=X(IMP1 ,J) 

DO 10 IS=1, NSEAS 
DO 10 IL=1,NLEV 
DO 10 IV=1,NVAR 

READ(8) ((VAR(LON,LAT,IV,IL,IS), LON=l,IMPl),LAT=l,JNP) 

C FOR ZONAL AVERAGES 

READ(8) (V ARZ(LAT,IV,IL,IS),LAT=1,JNP) 

10 CONTINUE 


XIX 





i. E-P FLUXES 


W3RWH.ECMWF.EPFLUX.SNMN.LOW 

Low pass (v'v'-u'u')/2 at 200mb, - u'v’ at 200mb, vT' a t850mb 
W3RWH.ECMWF.EPFLUX.SNMN.BAND 

Band pass ( v ' v ' “ u ' u ’) / 2 at 200mb, -u’v' at 200mb, V ' T ' at 850mb 
W3RWH.ECMWF.EPFLUX.SNMN.UNFILT 

Unfiltered ( vv - u ' u ') ^ 2 at 200mb, -u * v ' a t 200mb, V 'T' at 850mb 


READ FORMAT 

C NSEAS IS THE NUMBER OF SEASONS 
C IS=1 DJF (AVERAGED FOR 1980-87) 

C IS=2 MAM (AVERAGED FOR 1980-87) 

C IS=3 JJA (AVERAGED FOR 1980-87) 

C IS =4 SON (AVERAGED FOR 1980-87) 

C IS =5 ANNUAL AVERAGE (AVERAGED FOR 1980-87) 

C NVAR IS THE NUMBER OF VARIABLES IN THE PARTICULAR DATA SET 

PARAMETER(IMP1=73,JNP=46,NSEAS=5,NVAR=3) 

DIMENSION VAR(IMP1,JNP,NVAR,NSEAS) 

C J=1 IS THE SOUTH POLE, 1=1 IS THE DATELINE 
C FIELDS ARE WRAPPED SO THAT X(1,J)=X(IMP1,J) 

DO 10 IS=1, NSEAS 
DO 10 IV=1,NVAR 

READ(8) ((V A R(LON,LAT,I V,IS ) , LON=l,IMPl),LAT=l,JNP) 

10 CONTINUE 


xx 



ii. OTHER GLOBAL FIELDS 


W3RWH.ECMWF.VARIANCE.SNMN.LOW 

Low pass VZ'Z' at 500mb, (v'v’ + u'u') at200mb 

W3RWH.ECMWF.COVARNCE.SNMN.LOW 

Low pass u'v’ at 200mb, u'T' at 850mb, v'T' at 850mb, u’q’ a t 850mb, Vq’ a t 850mb 
W3RWH.ECMWF.VARIANCE.SNMN.BAND 
Bandpass VZ'Z' at500mb, (v'v'+u'u') at200mb 
W3RWH.ECMWF.COVARNCE.SNMN.BAND 

Band pass u'v' at 200mb, uT at 850mb, 7T at 850mb, u’q’ at 850mb, v 7 ^ at 850mb 

W3RWH.ECMWF.VARIANCE.SNMN.UNFILT 
Unfiltered Vz 7 ^ at500mb, (vV + uV) a t200mb 
W3RWH.ECMWF.COVARNCE.SNMN.UNFILT 

Unfiltered u'v' at 200mb, u'T' at 850mb, v'T' at 850mb, u'q' a t 850mb, v’q' a t 850mb 

READ FORMAT 

C NSEAS IS THE NUMBER OF SEASONS 
C IS=1 DJF (AVERAGED FOR 1980-87) 

C IS =2 MAM (AVERAGED FOR 1980-87) 

C IS=3 JJA (AVERAGED FOR 1980-87) 

C IS =4 SON (AVERAGED FOR 1980-87) 

C IS =5 ANNUAL AVERAGE (AVERAGED FOR 1980-87) 

C NVAR IS THE NUMBER OF VARIABLES IN THE PARTICULAR DATA SET 

PARAMETER(IMP1=73,JNP=46,NSEAS=5,NVAR=3) 

DIMENSION VAR(IMP1,JNP, NVAR, NSEAS) 

C J=1 IS THE SOUTH POLE, 1=1 IS THE DATELINE 
C FIELDS ARE WRAPPED SO THAT X(1,J)=X(IMP1,J) 

DO 10 IS=1, NSEAS 
DO 10 IV=1,NVAR 

READ(8) ((VAR(LON,LAT,IV,IS), LON=l,IMPl),LAT=l,JNP) 

10 CONTINUE 


XXI 



iii. ZONAL AVERAGES 


W3RWH.ECMWF.ZONAL.FLUX.SNMN.LOW 
Low pass [u'v 1 ], [vT] 

W3RWH.ECMWF.ZONAL.FLUX.SNMN.BAND 
Bandpass [u'v'], [ V 'T'] 

W3RWH.ECMWF.ZONAL.FLUX.SNMN.UNFILT 
Unfiltered [u'v 1 ], [v'T'] 

READ FORMAT 

C NSEAS IS THE NUMBER OF SEASONS 
C MOMENTUM FLUXES ARE AVERAGES FOR 1980-87 
C HEAT FLUXES ARE AVERAGES FOR 1980-85 
C IS=1 DJF 

C IS=2 MAM 

C IS=3 JJA 

C IS=4 SON 

C IS =5 ANNUAL AVERAGE 

C NVAR IS THE NUMBER OF VARIABLES IN THE PARTICULAR DATA SET 

PARAMETER(JNP=46,NLEV=7,NSEAS=5,NVAR=2) 

DIMENSION VAR(JNP,NVAR,NSEAS) 

C J=1 IS THE SOUTH POLE 

DO 10 IS=1, NSEAS 
DO 10 IV=1,NVAR 

READ(8) ((VAR(LAT,ILEV > IV,IS),LAT=l,JNP),ILEV=l t NLEV) 

10 CONTINUE 


xxii 



C) TIME SERIES OF ZQNALLY AVERAGE MONTHLY MEANS 


These fields are stored in a partitioned data set 


W3CKP.ECMWF. ATLAS (member) 


where the members are 


TVAR: Transient variance statistics (eg. u’u') 
ZVAR: Standing variance statistics (eg. u*u*) 
TCOV: Transient covariance statistics (eg. u'v') 
ZCOV: Standing covariance statistics (eg. u*v*) 


READ FORMAT 

PARAMETER(NMONS=96,JNP=46,NVAR=4) 
DIMENSION VAR(NMONSJNP) 


C NVAR IS THE NUMBER OF VARIABLES = 6 (VARIANCE) OR 4 
(COVARIANCE) 




VARIANCE 

COVARIANCE 

c 

IV=1 

(UU) 1/2 AT 200 MB 

UV AT 200 MB 

c 

IV =2 

(VV)‘« AT 200 MB 

VT AT 850 MB 

c 

IV =3 

(TT)i/2 AT 850 MB 

VQ AT 850 MB 

c 

IV=4 

(TT)'« AT 200 MB 

VZ AT 200 MB 

c 

rv=5 

(Q Q)>« AT 850 MB 


c 

IV =6 

(ZZ) l « AT 500 MB 



C J=1 IS THE SOUTH POLE, 1=1 IS JAN OF 1980, 1=96 IS DEC OF 1987 
C 

DO 10 IV=1,NVAR 
READ(8) IV, VAR 
10 CONTINUE 


xxm 






VII. THE SECOND MOMENT FIELDS 
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latitude latitude 
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